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ABSTRACT 

The pancake collapse of pregalactic clouds under UV background radiation 
is explored with a one-dimensional sheet model. Here, attention is concentrated 
on elucidating the basic physics on the thermal evolution of pregalactic clouds 
exposed to diffuse UV radiation. So, we treat accurately the radiation transfer 
for the ionizing photons, with solving chemical reactions regarding hydrogen 
molecules as well as atoms. The self-shielding against UV radiation by H 2 Lyman- 
Werner bands, which regulates the photo-dissociation of hydrogen molecules, is 
also taken into account. As a result, it is found that when the UV background 
radiation is at a level of 10 _22 (z//z/i) _1 erg s _1 cm _2 Hz _1 str _1 , the cloud evolution 
bifurcates with a critical mass as M SB = 2.2 x lO n M [(l + 2 c )/5] , where 
z c is the final collapse epoch. A cloud more massive than M$b cools below 
5 x 10 3 K due to H_2 line emission at the pancake collapse and would undergo 
the initial starburst. The pancake possibly evolves into a virialized system in 
a dissipationless fashion. Consequently, this leads to the dissipationless galaxy 
formation at 3 ^ z c ^ 10. A cloud less massive than Mg B cannot cool by 
H 2 emission shortly after the pancake collapse, but could cool in the course 
of shrinking to the rotation barrier. This is likely to lead to the dissipational 
galaxy formation at relatively low redshifts as ^ z c ^ 4. The present results 
provide a solid physical mechanism which controls the star formation efficiency 
in the pregalactic clouds. In the context of a standard CDM cosmology, M$b lies 
between la and 2a density fluctuations. 

Subject headings: galaxies: formation — radiation transfer — shock heating - 
molecular processes 
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1. INTRODUCTION 

It has been widely accepted that the formation of the first generation of objects, say 
Pop III objects, is regulated by the cooling by primordial hydrogen molecules. A number 
of authors have explored the formation of Pop III objects in so-called 'dark age' (z > 5) by 
concentrating on the role of hydrogen molecules (Matsuda, Sato, & Takeda 1965; Yoneyama 
1972; Hutchins 1976; Carlberg 1981; Palla, Salpeter, & Stahler 1983; Susa, Uehara, & Nishi 
1996; Uehara et al. 1996; Annonis & Norman 1996; Tegmark et al. 1997; Nakamura 
& Umemura 1999). Also, even if the gas contains metals, the metallic line cooling is over- 
whelmed by the hydrogen/helium cooling at T > 10 4 K when the metallicity is lower than 
10~ 2 Z Q (Bohringer & Hensler 1989). In fact, recently the metal abundance in the inter- 
galactic space is inferred to be at a level of 10 _3 Z Q (Cowie et al. 1995; Songaila & Cowie 
1996; Songaila 1997; Cowie & Songaila 1998). Moreover, at lower temperature of T < 10 4 K, 
the H 2 cooling is estimated to be still predominant as long as the metallicity is lower than 
10~ 2 Z Q , from the comparison of the cooling of the solar abundance gas (Spitzer 1978) with 
the maximal H 2 cooling (Kang & Shapiro 1992). Hence, hydrogen molecules are likely to 
play a fundamental role on the formation of galaxies under the metal poor environment. 

In addition, there seem to be situations where UV background radiation significantly 
influences the dynamical as well as thermal evolution of pregalactic clouds. At the epochs of 
z < 5, the existence of UV background of lCT 21±0 - 5 erg s _1 cm _2 str _1 Hz _1 at the Lyman limit 
is inferred from so-called proximity effect of Lya forest (Bajtlik, Duncan, & Ostriker 1988; 
Giallongo et al. 1996). Also, the Gunn- Peterson optical depths show that the intergalactic 
space was in fact highly ionized (Gunn & Peterson 1965; Schneider, Schmidt, & Gunn 1989, 
1991). The UV background is attributed to the radiation from quasars and partially from 
young galaxies (e.g. Giallongo et al. 1996). Even in the dark age of the universe, say z > 5, 
the first generation of objects as well as unseen proto-quasars might be external photoion- 
ization sources for subsequently collapsing objects. So far, many authors have considered 
the dynamical effects produced by the UV background radiation (Umemura & Ikeuchi 1984; 
Dekel & Rees 1987; Babul & Rees 1992; Efstathiou 1992; Chiba & Nath 1994; Thoul & 
Weinberg 1996; Quinn, Katz, & Efstathiou 1996; Navarro & Steinmetz 1997), to account 
for Lja clouds or to reconcile the paradox that in the hierarchical clustering paradigm for 
galaxy formation, low mass galaxies are overproduced compared with observations (White 
& Frenk 1991; Kauffman, White, & Guiderdoni 1993; Cole et al. 1994). The simulations 
have been hitherto based upon the assumption that the medium is optically thin against 
UV photons. However, the gas clouds become optically thick as the gravitational collapse 
proceeds. The significance of the optical depth has been stressed by the accurate treatment 
on radiative transfer (Tajiri & Umemura 1998). 
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The primary effects of UV radiation are photoionization, UV heating, and photo- 
dissociation of H 2 [Haiman, Rees, & Loeb 1996 (hereafter HRL), 1997; Corbelli, Galli, & 
Palla 1997; Kepner, Babul, & Spergel 1997]. HRL, with including radiation transfer effects 
in one-zone clouds, estimated the UV heating rate and the cooling rate under the assumption 
of the chemical equilibrium. They found that the cooling rate can exceed the heating rate at 
T ^ 10 4 K for high density. However, basically the chemical reactions regarding H 2 formation 
are not in equilibrium (e.g. Susa et al. 1998; Bertoldi & Draine 1996; Diaz-Miller, Franco & 
Shore 1998; Hollenbach & Tielens 1999). In particular, if ionization processes such as shock 
heating or photoionization take place, the chemical processes become out of equilibrium (e.g. 
Shapiro & Kang 1987). Thus, the non-equilibrium process for the formation of H 2 should be 
involved, coupled with hydrodynamical process. 

In this paper, we re-examine the thermal and dynamical evolution of pregalactic clouds 
under UV background radiation. In the present analysis, we solve properly the radiative 
transfer of diffuse UV photons and include non-equilibrium chemical reactions regarding 
H 2 formation due to shock-ionization and photoionization. The goal is to elucidate the 
key physics on the effects of the UV background radiation upon the thermal evolution of 
pregalactic clouds. For the purpose, we make the dynamical model as simple as possible. 
A cosmological density perturbation far beyond the Jeans scale forms a flat pancake-like 
disk. Although the pancake formation is originally studied by Zel'dovich (1970) in the 
context of the adiabatic fluctuations in baryon or hot dark matter-dominated universes, 
recent numerical simulations have shown that such pancake structures also emerge in a 
CDM cosmology (e.g., Cen et al. 1994). Thus, the pancakes are thought to be a ubiquitous 
feature in gravitational instability scenarios. Here, we consider pancakes in the plane-parallel 
symmetry. Also, previous 1D/2D/3D simulations on the pancake collapse in a dark matter- 
dominated universe show that the pancake is finally dominated by baryons after the caustics 
(Shapiro, Struck-Marcell, & Melott 1983; Bond et al.1984; Shapiro & Struck-Marcell 1985; 
Yuan, Centrella, & Norman 1991; Umemura 1993; Cen et al. 1994; Miralda-Escude & 
Rees 1994; Hernquist et al. 1995; Zhang, Anninos, & Norman 1995; Anninos & Norman 
1996). The presence of dark matter increases the shock velocity of the falling matter into 
the caustics. According to previous analyses (e.g. Shapiro & Kang 1987), such shocks 
change the thermal evolution in the absence of a UV background when the temperature 
exceeds 10 4 K since the temperature controls H 2 fraction in the postshock regions. However, 
under UV background radiation, the difference of shock velocity will not alter the results 
significantly, because the sheets are initially ionized by UV radiation and therefore quickly 
heated up to ~ 10 4 K. Thus, the thermal evolution of shocked region will be rather similar 
for any infall velocity, regardless of the shock temperature. Thus, in order to highlight the 
relevant physics of thermal processes, we deliberately exclude the dark matter contribution. 
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Nonetheless, some dynamical effects are anticipated from dark matter, e.g. Jeans instability. 
They are discussed later in the light of the present results. 

In section 2, we describe the basic equations and initial conditions. In section ^ the 
numerical results are presented. In section |], we make the physical interpretation of numer- 
ical results and present the condition on the cooling of a collapsing pancake. Section [5] is 
devoted to the implications for galaxy formation under UV background radiation. In the 
last section, we summarize the results. 



2. Formulation 

2.1. Basic Equations 

In this section, we give the basic equations for hydro dynamical calculations. We assume 
the plane-parallel symmetry throughout this paper. A set of the ordinary hydrodynamical 
equations is given as below. 

d (1/p) du = Q 
dt dm 

- + = U9 _ A( , r , B) , (3 ) 

where p,p,u, and z denote the density, the thermal pressure, the velocity, and the length 
measured from the mid-plane, respectively, and E is the total energy per unit mass, 

E - ^4 

m(z) is the column density, 

m(z) = / p(z)dz, (5) 
J o 

and g is the gravitational acceleration, 

g = 2nGm(z). (6) 

The symbol A(p, T, y^) denotes the cooling rate minus heating rate per unit mass due to the 
radiative and chemical processes of i-th species. We take H 2 cooling rate from Hollenbach & 
Mckee (1979), and the other atomic cooling rates are the same as those in Shapiro & Kang 
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(1987), except the cooling rate concerning on helium. Photoheating rate is calculated by 
taking into account the radiation transfer effects. The formulation is given in Appendix [A[ 
The equation of state is assumed as p = pkT/fiM, where k denote the Boltzman constant, 
and hm is the mean molecular weight of the five species discussed below in this section. The 
hydrodynamical equations are solved by the Piecewise Parabolic Method (PPM) described 
in Colella & Woodward (1984). PPM is one of the accurate methods to resolve the strong 
shock front. We take typically 800 spatial Lagrange grids. The hydrodynamical routine is 
tested by the Sod's problem. This code resolves the shock front typically by a few mesh. 

In the equations (fB), (0), and ([3]), the cosmic expansion is neglected. It is because we 
have considered the massive clouds far beyond the Jeans scale (Aj ~ 10 [(1 + z)/10] _1 ' 5 kpc) 
and also postulated that the initial stage is close to the maximum expansion for a density 
fluctuation. The effects of the cosmic expansion could be important for a density perturba- 
tion near the Jeans scale, since the temperature is raised up to ~ 10 4 K by UV background 
radiation in spite of the expansion and consequently the dynamics could be significantly 
influenced by the UV heating (e.g. Umemura & Ikeuchi 1984). However, far beyond the 
Jeans scale, the thermal pressure of gas around 10 4 K anyway does not affect the dynamical 
evolution until the maximum expansion as shown by previous calculations (e.g. Thoul & 
Weinberg 1996). Therefore, we attempt to pursue the history after the maximum expansion. 

The hydrodynamical equations are coupled with the non-equilibrium rate equations for 
chemical reactions: 

-JT = k jyj + n Yl hiVkVi + n 2 k mns y m y n y s , (7) 

°^ j k,l m,n,s 

where yi is the fraction of i-th species, yi = n«/n, with n being the number density of 
hydrogen nuclei, and fc's are the coefficients of reaction rates (see Table |l|). In the present 
calculations, we take into account six species e, H, H + , H~, H 2 , and H^. Here, we neglect 
helium. This is because the helium lines are not an effective coolant for lower temperature 
(T ^ 10 4 K) in which we are here especially interested. Nonetheless, some of helium lines 
including recombination lines may affect the ionization structure, because they are energetic 
enough to ionize HI. However, according to Osterborck (1989), the difference between the 
ionization structure of pure hydrogen and hydrogen + helium gas for 40000 K black body 
radiation is small. In the case of the power-law spectrum of UV background radiation which 
we consider, the difference could be larger, because the higher energy photons contribute 
more to the ionization of helium compared to the case of black body radiation. In fact, in 
the recent calculation by Abel & Haehnelt (1999), in which the radiation transfer effects of 
helium recombination photons are taken into account, the temperature of the clouds could 
differ by a factor of 2. Although we have dismissed helium to focus upon elucidating the 
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key physics, we should keep in mind that helium could be important in highly quantitative 
arguments. 

The photoionization and heating processes due to the external UV radiation are pursued 
by solving the frequency- dependent radiative transfer equation for hydrogen, 



where /i, k u , and i] v denote the direction cosine, the absorption coefficient, and the emissivity, 
respectively. The frequency-integration in the transfer equation can be done analytically for 
the continuum, while the recombination line near the Lyman limit is separately treated by 
solving transfer equation with a source term which comes from the recombination. The 
details are described in Appendix [A]. Equation (§) assumes steadiness of the radiation 
fields. When a neutral cloud is exposed to background ionizing radiation, ionized regions 
spread inward. In the ionized regions, the radiation propagates basically at the light speed, 
which is typically 10 3 times larger than the hydrodynamical velocity of the system. Thus, 
as for the radiation fields, it is sufficient to solve the steady radiative transfer equation. 
But, the ionization front propagates with a different speed, which is determined by the 
balance between the number of neutral atoms flowing through the front per second and 
the corresponding number of ionizing photons reaching the front. Resultantly, the speed of 
ionization front is much lower than the light speed, typically 500 km s _1 (ra/10 _2 cm~ 3 ) -1 
for UV background of 10 _21 erg s _1 cm _2 str~ 1 Hz~ 1 . In order to capture the ionization front 
propagation, we take the following numerical procedures. First, the steady radiative transfer 
equation (H) is solved, where typically 800 spatial grids are used, which number is equivalent 
to the number of Lagrange meshes assigned for hydrodynamics, and 40 angular grids are 
adopted. We have checked that the resultant accuracy is not significantly changed if we use 
the larger number of angular grids. The integration of equation ([8]) is performed with the 
summation of the exact formal solution in every grid (Stone, Miharas, & Norman 1992). This 
method is accurate, even if the optical depth of each grid can exceeds unity. Also, the validity 
of this method is enhanced by the fact that the ionization front is not sharply edged, but fairly 
dull for a power law of UV radiation (Tajiri & Umemura 1998). For the photons which cause 
H2 dissociation through the Solomon process (e.g. HRL1996), we employ the self-shielding 
function given by Draine & Bertoldi (1996), instead of directly solving the radiation transfer 
equation. Next, the rate equations including photoionization are solved implicitly for a time 
step which is determined as bellow. Finally, the hydrodynamical equations with the thermal 
equation is integrated explicitly. 

The time step of the calculation At is taken as follows: 




8 



At 



mm 



(0.99t hydro 



0.8t grav ,0.3t cool ,0.3r eat ). 



(9) 
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Here £ h y dro represents the hydrodynamical time, and 0.99t hydro is employed so that the 
Courant condition is satisfied. t grav denotes the time scale of gravity, which is defined as 
^grav = \fpfw/\g\, where p, 7, p, and g represents the pressure, the adiabatic index, the 
density, and the gravitational acceleration, respectively t c ° o1 and t heat represent local cool- 
ing time, and heating time scale, respectively. With these procedures, we have checked the 
propagation speed of the ionization front defined by the point ?/hi = 0.1. In all of our sim- 
ulations, the I-front propagation speeds are less than 1 percent of the speed of light. Thus, 
the assumption of the steady radiation transfer is valid. 



2.2. Initial Conditions 

We assume the initial density distributions to have a cosine profile; 

n- mi (z) = n- mi [1 + ecos (nz/X)] , (10) 

where n ini and A are the initial mean number density and the thickness of the sheet (See 
Fig.l), and e denotes the density contrast, which is set to be 0.5 throughout this paper. The 
ranges of n m i and A which we consider are respectively 10 _6 cm~ 3 < n m \ < 5 x 10~ 2 cm~ 3 
and 0.3kpc < A < 3Mpc. The initial velocity is null for every mass layer, because the 
initial stage is implicitly assumed to be close to the maximum expansion stage of a density 
fluctuation. In general, an overdense region at the maximum expansion is surrounded by 
an underdense region which is still expanding. Hence, if we introduce the initial velocity 
distributions more realistically, the further mass may accrete from the underdense region 
after the pancake collapse. However, the envelope does not seem to contribute to the self- 
shielding because of its low density, and thus the self-shielding is determined by the column 
density of the first collapsed sheet. Hence, we neglect the surrounding expanding underdense 
region, and model only the overdense region. Moreover, the velocity gradient may exist in a 
collapsing overdense region, and then leads to the delay of mass accretion from outer regions. 
The degree of self-shielding is determined by the total recombination number per unit time 
(./Vr ec ) in the volume against the UV photon number per unit time (iVuv) from the boundary. 
If the envelope has a density distribution in proportion to r x , N rec is proportional to r 2x+1 , 
whereas iVuv is constant, as far as the envelope undergoes nearly sheet-like collapse. Thus, an 
envelope with a distribution steeper than r -1 / 2 would not contribute to the self-shielding. In 
other words, the inner regions which collapse nearly simultaneously contribute mainly to the 
self-shielding. Although there are such qualitative expectations, we cannot assess accurately 
the effects by the initial velocity gradient and the cosmological expansion, because we have 
not included them. They should be quantitatively investigated in the future analyses. 

The initial temperature and chemical composition are given respectively by the thermal 
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equilibrium and by the ionization equilibrium. Also, the fractions of H 2 , H _ and are 
initially determined by the chemical equilibrium. When determining the initial equilibrium 
state, H2 cooling rate is ignored, because we set the initial condition so that the UV heating 
should overwhelm the H2 cooling. It will be shown below that this condition can be satisfied 
even if the cloud is almost neutral. 

The incident radiation intensity, PJ 1 , is assumed to have a power low spectrum as 

PJ 1 = I 21 (v/v L y a lO- 21 erg s^cm^str^Hz" 1 , (11) 

where vl is the frequency at the Lyman limit. The observations require the diffuse UV 
radiation to be at a level of J 2 i = 10 ±0 5 at z — 1.7 — 4.1 (Bajtlik, Duncan, & Ostriker 1988; 
Giallongo et al. 1996). At higher redshifts, we have no indication for J 2 i- Here, we assume 
I21 = 0.1, and the scaling to PJ 1 is argued. As for a, two typical cases are considered, i.e. 
a = 1 or 5, which are the same choice as Thoul & Weinberg (1996). The former represents 
a quasar-like spectrum, and the latter resembles a spectrum of star-forming galaxies around 
the Lyman limit frequency. 

3. Thermal and Dynamical Evolution 

In this section, we show the characteristic behaviors of the thermal/dynamical evolution 
resulting from numerical runs, and elucidate the dependence of the final states upon the 
initial states. 



3.1. Hard Spectrum Case: a — 1 

Figures 2 and 3 show the time evolution of spatial distributions of various physical 
quantities in the case of a — 1. Figure 2 represents the case with A = 100 kpc and n ini = 10~ 2 
cm -3 . In this figure, each horizontal axis represents the Lagrange coordinate, normalized by 
the half total surface density, m . The thick line shows the distribution at t — 1.28/ \/4nGpi ni , 
where p in \ = m p n in i with m p being the proton mass. The calculations are terminated there. 
Figure 2 shows that in the course of the collapse, a thin cold sheet quickly forms due to the 
efficient molecular cooling, which is confined by a shock-heated high temperature layer. The 
lower panels in Figure 2 show the corresponding HI and H2 distributions. Only an outer 
envelope (m/m 0.9) is photoionized and the rest parts of the sheet are self-shielded against 
UV radiation. As a result, the sheet is mostly neutral throughout the dynamical evolution, 
excepting the collisionally ionized postshock regions. In the early phase, UV heating exceeds 
the H 2 line cooling even in the almost neutral regions. The H 2 cooling becomes efficient as 
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the collapse proceeds to a further degree, and finally the temperature drops down to ~ 100 
K. Then, the collisional dissociation of hydrogen molecules becomes ineffective, and thus the 
fraction of hydrogen molecules at m = continues to increase. 

Figure 3 represents the case with lower density and larger size, i.e., n ini = 5 x 10~ 5 cm~ 3 
and A = 1 Mpc. The sheet also collapses because the initial thickness is larger than the 
initial Jeans length. However, the central temperature does not drop below ~ 10 4 K, and 
the resultant H 2 fraction is too low to overwhelm the UV heating. Why the evolution results 
in the different final states from the case shown in Figure 2 is basically understood by the 
difference of the degree of self-shielding against the external UV radiation. In the previous 
case, the most regions are so strongly shielded from the external UV that abundant H 2 form 
and the cooling rate due to H 2 becomes predominant. Contrastively, in the present case, 
the external UV permeates the inner regions to a considerable degree, so that UV heating 
overwhelms H 2 cooling everywhere until the collapse stops due to the thermal pressure. 



Figure 4 shows the case of a softer spectrum, a — 5. The other initial parameters are 
the same as Figure 3. Compared to the previous hard spectrum cases, the inner parts of the 
sheet are quickly self-shielded as the collapse proceeds. This is due to the fact that there are 
fewer high energy photons which penetrate the deep inside the sheet owing to the smaller 
ionization cross-section (see e.g. Tajiri & Umemura 1998). As a result, the photoionization 
rate as well as the UV heating rate damps strongly in the inner regions, so that a cold layer 
(T < 10 4 K) due to H 2 cooling emerges eventually. The self-shielded regions grow thicker 
and the the fraction of H 2 goes up to a level of 10~ 3 similar to the case of Figure 2. 



In Figure 5, the final physical states of the sheets for a — 1 at t — 1.28/y/AnGpi a i are 
summarized on the initial density-to-size diagram. Depending upon the initial thickness A 
and mean density rii n i, the sheets evolve in different ways. Filled circles in Figure 5 denote 
the initial conditions on which the clouds eventually cool below 5 x 10 3 K during the collapse. 
Open circles denote the ones on which the clouds cannot cool (T > 5 x 10 3 K). The boundary 
of the evolutionary bifurcation is well fitted by a simple power low relation as 



3.2. Soft Spectrum Case: a = 5 



3.3. Summary of the Dependence on Initial States 




(12) 
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The small open squares, open triangles, and filled triangles in Figure 5 trace the evolutionary 
paths for three typical parameters. In the paths, the time-dependent thickness (H) is defined 
by the thickness within which 80% of the total mass of the sheet is contained and the time- 
dependent mean density (n) is defined by the spatially averaged density over H; 



: / -dm, (13) 
Jo 



r 0.9m 1 

H = 2 

/o p 
2 r H 

n = — I ndz. (14) 



o 



The open squares and open triangles are respectively the results for the initial parameters 
of (A, fii n i) = (IMpc, 5 x 10~ 5 cm~ 3 ) and (lOOkpc, 1 x 10 _3 cm -3 ). For such initial conditions, 
the sheet shrinks due to Jeans instability [A > \j(T = 10 4 K)]. However, the collapse ceases 
by the thermal pressure before it intersects the critical line marked by A s hi e id, above which 
the sheet can be self-shielded and cool down due to formed H 2 . (The further details on A sh i e id 
are discussed in the next section.) On the other hand, the filled triangles, corresponding to 
the run with (A,n ini ) = (lOOkpc, 1 x 10~ 2 cm -3 ), come into the regions A > A s hi C id before the 
thermal pressure prevents the sheet from collapsing. Consequently, this sheet is shielded and 
cools below T ~ 5 x 10 3 K. 

Figure 6 shows the dependence on the initial states for a = 5. Because of the smaller 
number of runs, it is somewhat hard to divide clearly the parameter space. If we dare to 
draw the boundary, it is proportional to n^' 66 , which is numerically 



0.66 

A coo , = 6.0x 1 0-kpc( L0x ^ w3 ) (15) 



This criterion is slightly different from the previous case, in the sense that A coo i for a = 5 
does not cross the intersection point of Aj(T = 10 4 K) = A sh i C id, although A coo i for a = 1 
does. 



4. Physical Conditions for H 2 Cooling under UV background 

In this section, we attempt to understand the underlying physics of the present numerical 
results. To begin with, we introduce a scale, A s hi e id, which characterizes the relative efficiency 
of H 2 cooling to UV heating. Then, coupled with some physical arguments, we try to estimate 
the boundary of the bifurcation, A coo i, which has been derived by the present numerical 
calculations. 

As a measure of the penetration of UV photons into the cloud, a shielding length, A s hi e kb 
is defined by the balance between the UV heating rate and the H 2 cooling rate, which are 
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respectively 

Heating Rate = n(z)y m T m {t„ uv (z)) , (16) 
Cooling Rate = n(z) 2 A n2 {y H2 {z), T(z)), (17) 

where r H i is the UV heating rate per particle and A Ha denotes the H 2 cooling function. 
As shown in Appendix |A| |eq.( |A12| )l, in an optically thick regime, the UV heating rate is 
expressed by a power law of the optical depth at the Lyman limit, not by an exponential 
law. This is again because higher energy photons have the smaller ionization cross-section. 
The optical depth at the Lyman limit measured from the outer boundary (z = z ont ) is 



Zout 

dz ny m a. 



where a v is the ionization cross-section. If the slab is almost neutral, 



t Vl ,v(°) - / dzna UL , (19) 
Jo 

= n ini a UL X/2. (20) 



Equating fll6|) with fllTD at the midplane of the slab (z = 0), and also using equations (p0| 
and ( |A12| ), we can assess A s hi e id as 



n-m\Vv L \ 3n ini A H2 1 + pJ 

where (3 = 1 + (a — 1) /3. The A s hi e id is plotted in Figures 5 and 6 by long-dashed lines. The 
cooling rate has been assumed to be Ah 2 = 10 _26 cm 3 s _1 , which is a maximal value of H2 
cooling for the primordial gas that is once heated up to T ^ 10 4 K (Shapiro & Kang 1987; 
Susa et al. 1998). 

The obtained A s hi e id represents the degree of self-shielding against the external UV at 
the initial stages: (1) For A > A s hi e id(^ini), the sheet is promptly self-shielded enough to be 
cooled by the H 2 . (2) For A < Aghieid^im), the sheet is not self-shielded initially, so that the 
gas is heated up to T ^ 10 4 K due to predominant UV heating. A similar argument with 
Ashieid is also applicable to each dynamical stage. In the case, we should just replace the 
density n ini by n, and A by H, both of which are defined in section |3|. Then, we interpret 
Ashieid as the boundary beyond which the sheet is quickly self-shielded and consequently 
cools down owing to the efficient formation of H2. In fact, the clouds denoted by filled circles 
in Figure 5 satisfy the condition H > A s hi e id before the collapse is abruptly decelerated by 
approaching Aj(T = 10 4 K). Eventually, the clouds cool down below 5 x 10 3 K. On the 
other hand, the clouds denoted by open circles above Aj in Figure 5 result in warm sheets of 
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~ 10 4 K, because the clouds first meet Aj(T = 10 4 K) before intersect A sh i C id- If we assume 
that the sheet clouds evolve keeping the relation nH = const, as anticipated in the ideal 
sheet collapse, the boundary of the bifurcation in the space (n^, A) is expected to be a 
line which satisfies n^X = const, and meets the intersection point of A s hi C id = Aj. This 
prediction, A coo i oc n^, is close to the numerically obtained boundary, A coo i oc n^' s - The 
small difference of the dependence may come from the fact that the simulations contain the 
spatial structure of density and temperature, and therefore partial self-shielding. 



We remark that the Jeans length Aj is maximally scaled by yfie/^o i n the presence of 
dark matter, where is the baryon density parameter and f2 is the total density parameter. 
However, the baryonic component dominates the gravity of the sheets in the final phase of 
the collapse, although dark matter component does initially. Thus, what actually happens is 
the intermediate of these two extreme cases. Although calculations that include dark matter 
are necessary to evaluate the effect quantitatively, we leave them elsewhere. 



Based upon the present numerical results, we consider the context of the galaxy for- 
mation under UV background radiation. First of all, the cooling by atomic processes of 
primordial gas is essential for the formation of H 2 molecules which control the star forma- 
tion in primordial objects. The elaborate analyses by Rees & Ostriker (1977) and also by 
Blumenthal et al. (1984) show that the cooling mass is roughly constant almost regardless of 
the virial temperature and density, which is ~ 1O 12 M . If there is UV background radiation, 
the cooling mass could alter because the line cooling by H and He at ^ 10 5 K is seriously 
reduced (e.g. Thoul & Weinberg 1996). The mass, however, is basically determined by the 
cooling mechanisms at 10 5 K. At such temperature, the cooling in photoionized gas is 
dominated by collisional ionization, radiative recombination, and thermal bremsstrahlung, 
and possibly the Compton cooling at high redshifts. These mechanisms may potentially cool 
the clouds with ~ 1O 12 M down to several 10 4 K. In fact, it is shown by numerical calcula- 
tions that the clouds with ~ 1O 12 M can cool and collapse under UV background (Umemura 
& Ikeuchi 1984). Thus, we set here the upper mass of primordial galaxies to be 1O 12 M . 
Then, the evolution of pregalactic clouds under UV background radiation is discriminated 
into four categories; (a) promptly self-shielded clouds, (b) starburst pancakes, (c) retarded 
star-forming galaxies, and (d) expanding clouds. In Figure 7, the parameter regions corre- 
sponding to the four categories are shown for a — 1, with lines of equal mass, say, 1O 8 M 
and 1O 12 M . The upper abscissa is the collapse redshifts, z C) if the initial stage is assumed 
to be at the maximum expansion of a spherical top-hat density fluctuation. Q = 1, h = 0.5, 
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and fifth 2 = 0.02 are assumed in order to interpret the maximal expansion density into the 
collapse redshift z c , where h is the present Hubble constant in units of 100km s -1 Mpc -1 . In 
this figure, la, 2a, and 3a density fluctuations expected in a standard CDM are also shown. 
Some further details for each category are discussed in the following, restricting ourselves to 
the case a — 1. 



5.1. Promptly Shielded Clouds 

The evolution of promptly self-shielded clouds above A s hi e id is virtually equivalent to 
the evolution under no UV background radiation. The collapse redshifts of such clouds are 
expected to be z c ^ 10. The formation of primordial objects under no UV background has 
been hitherto extensively studied by numerous authors (Matsuda, Sato, & Takeda 1965; 
Yoneyama 1972; Hutchins 1976; Palla, Salpeter, & Stahler 1983; Susa, Uehara, & Nishi 
1996; Annonis & Norman 1996; Tegmark et al. 1997). As a result, we can expect the 
first generation of objects in the mass range of 1O 6 M at z c = 15 down toward 1O 3 M at 
z c > 200 (Tegmark et al. 1997). 

5.2. Initial Starbursts in Pancakes 

The region (b) in Figure 7 results in cold pancakes due to H_2 cooling, in which initial 
starbursts may take place. The instability of a shock-compressed layer has been discussed 
by several authors (Elmegreen & Elmegreen 1978; Vishniac 1983; Lubow Sz Pringle 1993; 
Whitworth et al. 1994; Yamada & Nishi 1998). According to Elmegreen & Elmegreen (1978), 
the fastest growing modes have the size of the sheet thickness and the fragmentation timescale 
is tf rag = (— 47rp(0)u> 2 ) -1 / 2 , where w 2 ~ —0.14 in the high external pressure limit. We 
estimate the line mass Zo at fragmentation by a condition tf rag = tdyn, where td yn = p(0)/p(0). 
Under an assumption that cold (T ^ 10 3 K) sheets fragment into filaments, the resultant 
line mass is tabulated in Table |2] with other properties, where / is the ratio of Iq to the 
critical line density l c = 2kT /m p G, and x s and x c denote the mass fraction of the shocked 
matter and cooled matter, respectively. It is noted that / is smaller than unity for the fastest 
growing modes. Hence, for the fragmentation, it is necessary for other growing modes to 
accumulate further mass so that the filaments would be super-critical. The super-critical 
filaments would eventually bear massive stars (Nakamura & Umemura 1999). 

We can estimate the minimum mass of the starburst pancakes by M SB = PimA 3 ool . Again, 
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if the initial stage is at the maximum expansion, then 



l + z c \~ 4 - 2 (Q B h 2 



) 



-1.4 



M SB 



2.2 x 10 n M 



5 J \ 0.02 



(22) 



The collapse epochs of the clouds in region (b) range from z c ~ 3 to z c ~ 10. The clouds may 
undergo further shrinking after the sheet-like collapse, because the rotation barrier is smaller 
by Ag pin than the maximum expansion size, where A sp i n is the dimensionless spin parameter 
which is peaked around 0.05 (Heavens & Peacock 1988). Consequently, the violent relaxation 
would take place in a collisionless fashion. Hence, the starburst pancakes highly possibly lead 
to the dissipationless galaxy formation. It is worthy noting that Msb lies between la and 
2cr fluctuations in a standard CDM scenario. Thus, they could form a relatively clustered 
population of galaxies. 



The clouds with the smaller mass than MgB cannot cool below 10 4 K at the first pancake 
collapse, but they might pass A s hi e id in the course of shrinking down to the rotation barrier, so 
that the clouds are self-shielded against external UV radiation to cool down by H2 molecules. 
Assuming A spin = 0.05, the marginal size, A rot , above which the clouds can pass A s hi e id is much 
smaller than size of the Jeans unstable clouds for z c > 0. Thus, no cloud can be stopped by 
the rotation barrier before the cloud is self-shielded. In the clouds in the region (c), the star 
formation is likely to be retarded until the pancake disk shrinks considerably. Therefore, the 
region (c) tends to result in the so-called dissipational galaxy formation (e.g. Larson 1976; 
Carlberg 1985; Katz & Gunn 1991). The collapse epochs are relatively later (0 ^ z c < 4), 
compared to the region (b), and the region (c) is corresponding to CDM fluctuations lower 
than la. Thus, it is likely that the region (c) leads to a less clustered population of galaxies. 
Futhermore, for low redshifts (z ^ 1), the intensity of UV background is no longer at a 
level of I21 = 0.1, but could be smaller by two orders of magnitude (Maloney 1993; Henry 
& Murthy 1993; Dove & Shull 1994). Since A s hi C id oc / 2 i for a = 1 as seen in (^), the 
clouds with z ^ 1 in the region (c) might cool quickly in the course of shrinking down to the 
rotation barrier. 



5.3. Retarded Galaxy Formation 



5.4. Expanding Clouds 



In the region (d), the clouds are no longer gravitationally bound because of the en- 
hancement of thermal pressure due to UV heating. They are relatively low mass systems. 
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They could be Lyman alpha absorption systems which are seen in QSO spectra (Umemura 
& Ikeuchi 1984, 1985; Bond, Szalay, & Silk 1988). 

6. CONCLUSIONS 

We have numerically explored the thermal and dynamical evolution of pregalactic clouds 
under UV background radiation. The plane-parallel collapse of primordial gas clouds has 
been pursued, including chemical reactions with respect to hydrogen molecules as well as 
atomic hydrogen. Also, the radiation transfer for the ionizing photons has been properly 
treated. As a result, it is found that the cloud evolution under UV background branches off 
into four categories in the initial parameter space of density and size. They are (a) promptly 
self-shielded clouds, which evolve into dense objects with ^ 1O 6 M at collapse redshifts 
(z c ) greater than 10, (b) starburst pancake clouds with the mass higher than M SB = 2.2 x 
1O 11 M [(1 + z c )/5]~ 4 ' 2 , which lead to the dissipationless galaxy formation at 3 ^ z c <> 10, (c) 
retarded star-forming galaxies with the mass lower than Msb, which undergo star formation 
in the course of shrinking down to the rotation barrier at ^ z c ^ 4, consequently leading 
to the dissipational galaxy formation, and (d) less massive expanding clouds, which could 
be detected as Lyman alpha absorbers in QSO spectra. If we assume a standard CDM 
cosmology, density fluctuations of \o-1o coincides with Msb- That is, fluctuations higher 
than 2a result in the dissipationless formation of massive galaxies at 3 ^ z c ^ 10, and 
eventually constitute a clustered population. From a further realistic point of view, there 
must be local subgalactic clumps growing in a pregalactic cloud in a bottom-up scenario 
like a CDM cosmology. They are likely to be self-shielded earlier than uniform components. 
Thus, M S b can be considered as a measure of determining the bulge-disk ratios (B/D) of 
formed galaxies. Above Msb, the B/D ratio should be larger, while the ratio becomes smaller 
below M SB . 
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the Ministry of Education, Science, Culture, and Sports, 09874055 (MU). 
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A. Frequency Integration for Photoionization and UV Heating Rate 

In this appendix, a method to integrate the frequency dependence in radiation transfer 
equation is presented (see also Tajiri & Umemura 1998). This method is applied to plane- 
parallel calculations in this paper. However, it is also potent for the 3D calculations on 
cosmic ionization problem (Nakamoto, Umemura & Susa 1999). The frequency-dependent 
radiation transfer equation is given by a general form as, 

n • V l v = -K v l v + r]u (Al) 

where n is the unit directional vector, k v and r\ v denote the opacity and the emissivity, 
respectively. For the gas composed of pure hydrogen, the emissivity is almost null for v 
ul + Aut, where Ai^t denotes the width of the Lyman limit emission resulting from radiative 
recombination; 

Au T = kT/hp, (A2) 

where T, k and hp denote the gas temperature, and Planck constant, respectively. Typically, 
Avt/vl ~ 0.1 for T = 10 4 K. As far as higher energy photons of v ^ z/^ + Aujn are concerned, 
the solution of equation ( |A1| ) is simply 

I v = I?exp(-r v (w)), (A3) 

where J™ denotes the incident intensity at the boundary. The optical depth t v (oj) in the 
solid angle uj is rewritten in terms of the optical depth at Lyman limit, t Vl {oj), 

t u (u)=t„ l (u)(v l /v) 3 , (A4) 

because the photoionization cross-section is a v ~ <j Ul (vl/v) 3 - Using the equations ( |A3|) and 
(|A3|), we obtain the photoionization rate coefficient for this energy range of photons, k^f, in 
terms of an integration as 

L 



k m = I I -r^Vududv, 

Ju L +Av T J nv 



h J 3r^(cu) 1+Q / 3 



l(l + a/3,T VL (u))du, (A5) 



« <7 IA K,.»l)> (A6) 

where we have assumed PJ 1 = I l „ L {v / 'vl)~°- and j(a, b) denotes the incomplete gamma func- 
tion. r VL)V denotes the optical depth measured in the vertical directions of a slab, which is 
already introduced in section |. It is noted that the integration in frequency space has been 
carried out just analytically with the incomplete gamma function. In other words, we do 
not have to solve the radiation transfer equation for v ^ + Avt- 
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On the other hand, photons with Ul < v < vl + Av are scattered to produce diffuse radi- 
ation. For hydrogen, the scattering albedo is given by [«a(T) — a B (T)]/a J 4(T), where cxa{T) 
is the total recombination coefficient to all bound levels and a#(T) is the recombination 
coefficient to all excited levels. Hence, we have to integrate the radiation transfer equation 
including the emissivity by scatterings. In this frequency range, the frequency- dependence 
of the opacity and the emissivity are thought to be small, because Ai>t is ten times smaller 
than vi for T = 10 4 K. So, we approximate the opacity and emissivity to be constant in 
this frequency range. After the transfer equation is solved, we calculate the angle averaged 
intensity J UL in terms of the numerical solution I VL . Using J UL , we obtain the photoionization 
rate coefficient for this frequency range, fogf; 

k sca _ 47r Z^£Az/ T (A7) 



Finally, by adding ( |A6|) to (|A7| ) , we obtain the total ionization rate /cm; 



kg? + kg? . (A* 



The UV heating rate Thi is also obtained in a similar way. The result is, 

r H i = r^ + rgf, (A9) 



Tg{ = / JLfj^hiy - hu L )dudu, 

r7(l + (a - 1) /3, t Vl (u))hu L duj 



I in a 



h J 3r^(a;) 1 +( a - 1 )/ 3 
-hu L Hg?, (A10) 
T^T ~ ^J VL a VL Av T (All) 

In particular, for an infinite sheet, we have an asymptotic expression for the absorption 
part (rgf) in an optically thick limit as 

T a bs „ 2<^^r(l + (q-l)/3) ^.^3 
HI 3 2+(a-l)/3 " L ' v ' 1 ; 

where T(a) represents the gamma function. It should be noted that the UV heating rate 
is not proportional to exp(— t VL;V ), but is proportional to t^ 1 ^*-" -1 ^ 3 for large t VL;V . If we 
choose a = 1 as in the present paper, the UV heating rate is proportional to t~^ v . This 
shallow dependence upon the optical depth comes from a steep dependence of ionization 
cross section upon frequencies. 
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W ini 




Fig. 1. — The profile for the initial density distributions. The horizontal axis, z, is the height 
form the midplane. The vertical axis is the number density, n ini . n ini and A denote the mean 
number density of the slab and the wave length of the slab, respectively, e is the parameter 
which characterizes the initial density contrast (see the text). 
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Fig. 2. — The time evolution of the spatial distributions for various physical quantities is 
presented in the case of n m \ = 10~ 2 cm~ 3 , A = 100 kpc, and a — 1. The horizontal axis 
in each panel is the Lagrange coordinate normalized by the half total surface density, m . 
Note that m/mo = corresponds to the midplane of the sheet and m/rriQ = 1 is the surface. 
The upper left panel shows the density distributions and the upper right panel shows the 
temperature distributions. The lower left is the distributions for HI fraction, and the lower 
right is the distributions for H 2 fraction. The thick solid lines show the final stage, i.e., 
t = 1.28 x (AjcGpini)^ 1 ^ 2 , where p ini = m p n mi . 




Fig. 3. — Same as Figure 2, except that n- mi = 5 x 10 5 cm 3 , A = 1 Mpc. 




Fig. 4. — Same as Figure 2, except that a — 5. 
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io- 5 10- 4 10- 3 10- 2 0.1 

n ini , n cm- 3 

Fig. 5. — The physical states at t — 1.28/ \fAnGpi n \ with given initial parameters are summa- 
rized. The vertical axis and the horizontal axis denote the thickness and the mean density of 
the sheets. The open and filled circles represent the initial conditions for the runs which we 
performed. The open circles denote the sheets whose central temperature is high (T > 5 x 10 3 
K) after one dynamical time. The sheets denoted by the filled circles have cooled below 5x 10 3 
K. Solid line (A coo i) denotes the boundary of the cold sheets and the hot ones. Long dashed 
line (A s hieid) denotes the critical thickness above which the sheets are initially self-shielded 
against external UV heating. Short-dashed line is the Jeans length with T = 10 4 K. Small 
open squares, open triangles, and filled triangles represent the evolutionary sequences for 
three initial parameters (see text). 
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Fig. 6. — Same as Figure 5, except a = 5. 
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Fig. 7. — The initial parameter space is divided according to the evolution of clouds. Here, 
a — 1 is assumed. Solid line, long-dashed line, dot-dashed line, and short-dashed line are the 
same as those in Figure 5. The upper abscissa is the estimated redshifts of the final collapse, 
z c . The dot-dashed line can be translated into M SB = 2.2 x lO n M [(1 + z c )/5]" 4 - 2 , if the 
initial stage is assumed to be at the maximum expansion. Thin dot-dashed lines denote 
the constant mass corresponding to M — 10 8 and 1O 12 M . Assuming a standard CDM 
cosmology, la, 2a, and 3a density fluctuations are also shown. The evolution of pregalactic 
clouds under UV background radiation branches off into four categories; (a) first generation 
(promptly self-shielded clouds), (b) starburst (initial starburst pancakes), (c) retarded GF 
(retarded star formation in primordial galaxies), and (d) expanding clouds (unbound UV- 
heated clouds). 
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Table 1. Reaction rates 



Reactions References 

H++ e ^H + 7 Spitzer 1956 

H + e ^H"+ 7 de Jong 1972 

H~+ H ^H 2 +e Beiniek 1980 

3H ^H 2 + H Palla, Salpeter & Stahler 1983 

H 2 + H ^3H SK 

2H + H 2 ^2H 2 Palla, Salpeter & Stahler 1983 

2H 2 ^2H + H 2 SK 

H + c ^H++2e Lotz 1968 

2H ^H + H++ e Palla, Salpeter & Stahler 1983 

H + H+->H£+ 7 Ramaker & Peek 1976 

H ^H 2 + H+ Karpas, Anicich & Huntress 1979 

H 2 + H+->H£+ H Prasad & Huntress 1980 

H"+ H+^H 2 f + e Poularert et al. 1978 

Hj+ e ^2H Mitchell & Deveau 1983 

Hj+ H-^H + H 2 Prasad & Huntress 1980 

H"+ e ^H + 2e Duley 1984 

H _ + H ^2H + e Izotov & Kolensnik 1984 

H"+ H+^2H Duley 1984 

H + 7 ^H + + e see text 

H 2 + 7 ^2H Draine & Bertoldi 1996 



Table 2. Initial physical parameters and the properties of cylindrical fragments 



A kpc 


n ini cm" 


-3 


l Mepc- 1 


/ 


x s (shocked) 


x c (cooled) 


100 


1 x 10" 


-3 


8.7 x 10 3 


2.4 x lO^ 1 


0.54 


not cooled 


100 


2 x 10" 


-3 


6.3 x 10 2 


1.7 x 10" 1 


0.78 


0.47 


100 


1 x 10- 


-2 


1.5 x 10 1 


2.2 x 10~ 2 


0.36 


0.31 



